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I Abstract 
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The nucleon is modeled, using light front dynamics, as a relativistic system of three bound con- 
K*" i stituent quarks emersed in a cloud of pions. The pionic cloud is important for understanding low- 

momentum transfer physics, especially the neutron charge radius, but the quarks are dominant at high 
■ values of . The model achieves a very good description of existing data for the four electromagnetic 
, elastic form factors. 
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The recent exciting experimental results for the ratio proton elastic form factors Ge/Gm (or 
QF2/-^i)llli m and the impending high accuracy data for the neutron electric|3| and magnetic|§] form 
^ ' factors have re-ignited interest in the venerable goal of understanding the structure of the nucleon. 
^ ■ The aim of the present paper is to present a reasonable, workable model which describes currently 

available information and makes predictions testable against data taken at higher values of Q , or taken 
for improved accuracy. The model should have enough content so that its ultimate disagreement with 



X 



■ experiment elucidates some missing piece of physics. Poincare invariance and chiral symmetry are the 



principal tools used to construct the model. 

Poincare invariance is maintained or approximated by using light- front dynamics, in which fields 
are quantized at a fixed "time"=r = + = The r-development operator is then given 

by — = P~ . The canonical spatial variable is x~ = x^ — x^ , with a canonical momentum 
P+ = P^ + P^. The other coordinates are x_|_ and P^. The relation between energy and momentum 

of a free particle is given by: p~ = — , a relativistic kinetic energy which does not contain a square 
root operator. This allows the separation of center of mass and relative coordinates, so that the 
computed wave functions are frame independent. The use of the light front is particularly relevant for 
calculating form factors, which are probability amplitudes for an nucleon to absorb a four momentum 
q and remain a nucleon. The initial and final nucleons have different total momenta. This means that 
the final nucleon is boosted relative to the initial one, and therefore has a different wave function. The 
light front technique allows one to use boosts that are independent of interactions. 

We are concerned with the Dirac Fi and Pauli F2 (-^2(0) = k, the anomalous magnetic moment) 
form factors. The Sachs form factors are Ge = Fi — Q^/4M^F2, Gm = -^1 + ^2- We use the current's 
"good" component, J+, so that Fi{Q^) = (A^, j |J+| A^,!), QF2{Q^) = {-2Mn){N,] |J+| A^,i), with 
nucleon light-cone spinors, and in a frame with = and = = q^. 
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The model nucleoli consists of three relativistically moving, bound constituent quarks, which are 
surrounded by a cloud of pions. The quark aspects ||5|-||9|, will be discussed first. The original con- 
struction of this three-quark model was based on symmetry principles [^,||6|. The wave function is 
anti-symmetric, a function of relative momenta, independent of reference frame, and an eigenstate of 
the spin operator. Schlumpf[^ applied it to compute a variety of baryonic properties. Prank, Jennings 
& I used this model to predict a very strong decrease of Ge/Gm as a function of Q^, which has now 
been measured. Explaining the meaning of this result was left for a second paper in which imposing 
Poincare invariance was shown to lead to an analytic result that the ratio QF2 / Fi is constant for large 

and to a violation ||To[| of the helicity conservation rule. 

The wave function we use is given by 



(1) 



where ip is a spin-isospin color amplitude factor, the pi are expressed in terms of relative coordinates, 
the u{pi) are ordinary Dirac spinors, $ is a spatial wave function and the repeated indices pi are 
summed over. The specific form of ip is given in Eq. (12) of Ref. and earlier in Ref. The notation 
is that Pi = {pf,pi±). The total momentum is P = pi + p2 + P3. The relative coordinates are 
^ = p+/p+ r] = {pf +pt)/P'^, and k_L = (1 - Opi± - ^P2±, K_l = (1 - r/)(pi_L + P2±) - ??P3±- 

In computing a form factor, we take quark 3 to be the one struck by the photon. The value of 1 — 77 is 
not changed {q~^ = 0), so only one relative momentum, K_|_ is changed: K'_|_ = K_|_ — ?7q_L. We take the 
form of the spatial wave function from Schlumpf|^: <^(Afo 
operator for a non-interacting system: 
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, with Mq is the mass-squared 



M2 



+ 



+ 



m 



77(1 - 7?) r/^(l - 1-7? 



(2) 



Schlumpf's parameters are /? = 0.607 GeV, 7 = 3.5, m = 0.267 GeV. The value of 7 was chosen 
that Q^Gj^jiQ"^) is approximately constant for > 4 GeV^ in accord with experimental data. The 
parameter P helps govern the values of the perp-momenta allowed by the wave function $ and is 
closely related to the rms charge radius, and m is mainly determined by the magnetic moment of the 
proton. We shall use different values when including the pion cloud. 

The calculation of form factors is simplified by using completeness to express the wave function 
in terms of light cone spinors ul{p^ ,p, X), which are related to Dirac spinors by a unitary Melosh 
rotation evaluated in terms of Pauli spinors: \Xi), with (Aj|-Rj^J(p^)|s^) = UL{pi, Xi)u{pi, Si). Thus 
the wave function depends on Melosh rotated Pauli spinors: 



T Pi 



m + (1 - 7])Mo + i<T ■ (n X pa) 

(m + (l-77)A/o)2+p|^ 



(3) 



where the quantity in brackets is R\j{ps)- The spin-isospin wave function can then be thought of 
as constructed from the non-relativistic quark model, but with the replacement of Pauli spinors by 
those of Eq. (^. An important effect resides in the term (n x P3) which originates from the lower 
components of the Dirac spinors: the orbital angular momentum 7^ 0[pl|]. The term (n x P3) is 
also responsible for the flatness of the ratio F2{Q'^)/ Fi{Q'^). 

We turn now to neutron properties. The three-quark model for the proton respects charge sym- 
metry, invariance under the interchange of u and d quarks, so it contains a prediction, shown in Fig. |l| 
(compared with data from Ref. ||l^) for neutron form factors. We note that Geu would vanish in 
the non-relativistic limit, Rm 1, so the deviations from are solely due to relativistic effects. The 
resulting electric form factor, shown in the curve labeled relativistic quarks, is very small at low values 
of Q^, but at larger values of the prediction is larger than that of the Galster parametrization|13[| . 

The slope of Geu is related to the charge radius as GEniQ'^) — > —Q'^R^/^ with a measured value |14j| 
of i?^ = —0.113 lb 0.005 fm^. The three-quark model value is —0.025 fm^. To understand this small 
magnitude we express in terms of Fi^2 for small values of Q^. Then = R^ + Rp, where the Foldy 
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Figure 2: Diagrams 

contribution, R'p = 6Kn/4M^ = —0.111 fm^ is, by itself, in good agreement with the experimental 
data. But this does not guarantee success in explaining the charge radius because one needs to include 
the dependence of Fi which gives -Rf . In the three-quark model = +0.086 fm^ which nearly 
cancels the effects of Rp. Such a cancellation is a natural consequence of including the relativistic 
effects of the lower components of the quark Dirac spinors|15]. Another effect is needed. 



Sometimes a physical nucleon can be a bare nucleon emersed in a pion cloud. An incident photon 
can interact electromagnetically with a bare nucleon. Fig. with a nucleon while a pion is present. 
Fig. or with a charged pion in flight. Fig. |^. These effects are especially pronounced for the 
neutron G£;[^], at small values of Q^, because the quark effects are small. The tail of the negatively 
charged pion distribution extends far out into space, causing i?^ to be negative. Such contributions 



were computed long ago using the cloudy bag model[|l6[, which employed static nucleons. 

It is necessary to compute the effects of the pion cloud in a relativistic manner, to confront data 
taken at large Q^. This involves evaluating the Feynman diagrams of Fig. § using photon- nucleon 
form factors from our relativistic model, and using a relativistic vr-nucleon form factor. We define the 
resulting model as the light-front cloudy bag model LFCBM. The light-front treatment is implemented 
by doing the integral over the virtual pion four- momentum k^,'k± performing the integral over k~ 
analytically, re-expressing the remaining integrals in terms of relative variables (a = k~^/P~^), and 
shifting the relative j_ variable to L_l to simplify the numerators. Thus the Feynman graphs, Fig. 0, 



are represented by a single r-ordered diagram. The use of and the Yan identity|17] Sf{p) = 
J2s ^(P' s)u{p, s) I [j? — rr? + it) + ^'^ jlp^ allows one see that the nucleon current operators appearing 
in Figj2|b act between on-mass-shell spinors. 
The results can be stated as 

Fia(Q') = Z \Ff^{Q') + F^UQ") + Fi,,(g2)l , (4) 
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Table 1: Different parameter sets, units in terms of fm 



bet( legend) 


m A 7 -fin Hp 


1 solid 

2 dot-dash 

3 dash 


1.8 3.65 3.1 4.1 0.111 1.73 2.88 
1.7 3.4 3.1 3.9 0.110 1.79 2.95 
1.7 2.65 3.1 3.7 0.109 1.79 2.95 



where i = 1,2 denotes the Dirac and Pauli form factors, a = n,p determines the identity of the 
nucleon, and {Q'^) are the form factors computed in the absence of pionic effects. The wave 
function renormalization constant Z is determined from the condition the charge of the proton be 
unity: -Fip(Q^ = 0) = 1. Calculating the graph Fig. gives 

^ Jn 

+ (2F(°)(Q2) + i.(0)(Q2)) f («2q2/2)^ (5) 

^ Jn 

F2UQ') = {fSHq') + Fi^\Q')/2) j^m^o?) + + i^i?(Q')/2) /^(4a'M2 + 2i?). (6) 

where the integration measure is given by 

9' f d'L^^R{L^^^ \a)R{L'--^ \a), (7) 



N 



2(27r)3 J ^ a 



g is the vrN coupling constant, g /47r = 14 ,L]_ ' = L_Libaq_L/2, aD{k-^,a) = M a + kj_ + fi {1 — a), 
and R{k\,Q) = -j^^^"^- The vrN form factor is taken as 

F,;v(A;i,a) = exp[-(Z)(A;i,a)/2(l -a)A2)], (8) 



as used by Refs. [18| and [|^, and satisfies the constraints needed to maintain charge conservation po[| . 
Including the form factor this way uses the assumption that the form factor is an analytic function of 
k~ . The results (^J^ show that each term in the nucleon current operator contributes to both Fi and 
F2 ■ The evaluation of graph ^ yields 

FianiQ^) = -2F.(Q') / ^^^—^{aHl^ + Li - (1 - a)QV4) (9) 

J N Ot 

F2cn{Q^) = -2F^(Q2) ( ii^(2m2a(l - a)). (10) 

JN OL 

The proton form factors can be obtained by simply making the replacements n — > p in Eqs. (^,^ and 
-2 ^ +2 in Eqs. 

Eqs. ([4[-|l0|) completely specify the form of the calculation. But the LFCBM requires four param- 
eters m, /?, 7, A. Including pionic effects while continuing to use the original values of m, /3, 7, would 
lead to a satisfactory description of G£;„[|lT|, but would cause other computed observables to disagree 
with experiment. Thus a new set of parameters is needed. The following set of requirements is used to 
restrict the parameters. First, the magnetic moments of the proton and neutron must agree with mea- 
sured ones within 10%. We also require that the computed values of GMn(0.5), G£;„(l, 1.5), GMn(4), 
A*G'£;p/Gmp(5.5), Gmp(5.5) and Gmp(IO) agree with the measured values well enough so that the av- 
erage disagreement is about one error bar. There are many parameter sets that satisfy this criterion. 
We show results for three in Table 1[^ and in the figures. 

The first application of the LFCBM is to Geu and the results of using the three parameters sets of 
Table 1 are shown in Fig. ||. It is easy to find many parameters which provide a large pionic effect at 



4 



2 4 6 

Q' (GeV') 



10 



Figure 3: Gunl ^-Gd for the neutron. Data are from Ref. p2 



small values of . The agreement with existing data is good, and more higher quality data at larger 
values of is expectedQ. The next step is to compute Gun, which is expressed as Gmu/ fJ-Gn, where 
fi is the computed neutron mag netic moment and Gd = (l + /0.7lGeV'^)-'^ . The resuhs are shown 
in Fig. ^. The agreement between the present theory and existing data [^2| is excellent, but this will 
soon be tested by a new experiment Q. 

We turn now to the calculation of proton observables. Fig. |^ shows that the measured ratio of Dirac 
to Pauli form factors is reasonably well-reproduced. This is very similar to our earlier results ^, 
showing that the pion cloud effects are not very important for this ratio at relatively large values of 
Q^. These ratios are insensitive to the parameter set, and the results for sets 2 and 3 overlap. Thus, 
as stressed elsewhere ||9|, |ll| respecting symmetries is more important than including detailed dynamics 
in obtaining a constant ratio. Finally, the proton magnetic form factor, expressed as GM/fJ-Go, where 
fi is the computed proton magnetic moment, is shown in Fig. ^. For this case, set 3 seems to provide a 
"best" description of the present data[Q up to about = 20 GeV^. For higher values the calculation 
falls a bit below the data, perhaps indicating the need for the effects of perturbative QCD. 

These calculations show that the combination of Poincare invariance and chiral symmetry, as 
realized in the present light front cloudy bag model, is sufficient to describe the existing experimental 
data up to about = 20 GeV^. This is somewhat surprising as the model keeps onlytwo necessary 



effects. The effects of mixing of quark-configurations [ 24 1, the variation of the quark mass with (5^[|2| 
exchange currents [|2^ and an intermediate A|[l^ could enter, but these seem to have little infiuence 
within our model. 

Perhaps the strongest feature of the model is that it is testable in upcoming experiments. For the 
proton, QF2/F1 is predicted to be constant for values of up to about 20 GeV^. The neutron Geu, 
soon to be measured, is predicted as is its magnetic form factor also soon to be measured. 

There also are implications for other reactions. Hadron helicity conservation predicts that Q'^F2/Fi 
is constant [26 1 . This is not respected in present data, so there is no need to expect it to hold for a 
variety of exclusive reactions occurring at high < 5.5 GeV^. Examples include the large spin 
effects observed in pp elastic scattering ||28|] and the reaction — > np|29|, but there are many other 
possibilities. 
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